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Abstract 

The two long-baseline experiments T2K and NOvA have recently presented new findings. T2K has shown the first 
Ve appearance data while NOvA has released the first Vg appearance results. These data are of particular importance 
because they allow us to probe for the first time in a direct (or manifest) way the leptonic CP-violation. In fact, it is 
the first time that a hint of CP-violation arises from the comparison of the observations of neutrinos and antineutrinos. 
We consider the implications of such new results both for the standard 3-flavor framework and for the non-standard 
3-1-1 scheme involving one sterile neutrino species. The 3-flavor analysis shows a consolidation of the previous trends, 
namely a slight preference for sin <5 < 0, disfavoring CP conservation {6 - 0, tt) with a statistical significance close 
to 90% C.L., and a mild preference (at more than 68% C.L.) for the normal hierarchy. In a 3-1-1 framework, the data 
constrain two CP-phases ((5i3 s 6 and ^14), which exhibit a slight preference for the common value ^13 ^ ^14 ^ -njl. 
Interestingly, in the enlarged four neutrino scheme the preference for the normal hierarchy found within the 3-flavor 
framework completely disappears. This indicates that light sterile neutrinos may constitute a potential source of 
fragility in the capability of the two LBL experiments of discriminating the neutrino mass hierarchy. 


1. Introduction 


Neutrino physics is entering a new era. After the discovery of a relatively large value of @13, the determination of 
leptonic CP violation (CPV) and the identification of the neutrino mass hierarchy (NMH) have become realistic targets. 
In this respect, the long-baseline (LBL) experiments represent the first setups able to probe these properties. One of 
them (T2K), already operative since 2011, has provided the first observation of the conversion of muon neutrinos 
into electron neutrinos 111]. This result, in combination with the precise determination of the third mixing angle 0i3 


achieved with the dedicated reactor experiments (2, 


m, has provided the first indirect information on CPV. In fact. 
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all the recent global yl I 3 17|] and partial || 8 |, 19|] (including only reactor and T2K data) fits indicate a slight preference 
for CPV, with a best fit value of the fundamental CP-phase 6 close to the value -njl. 


A few weeks ago the T2K experiment has presented the first appearance results in the antineutrino channel J 10, 
1 in . In the same days the experiment NOvA has released the first appearance data in the neutrino channel lll2n . 


Although both measurements have a limited statistical power it is nonetheless interesting to consider their impact on 
the preexisting scenario. With this purpose, we here present both a 3-flavor and a 4-flavor analysis of the results of the 
two long-baseline experiments, updating our previous work ||8t]. 

We underline that the latest T2K and NOvA data considered in the present work are of particular importance. 
Arguably, they represent the most important piece of information in the current neutrino landscape. In fact, these data 
allow us to probe for the first time in a direct (or manifest) way the leptonic CPV. It is the first time that a hint of 
CPV arises from the comparison of the observations of neutrinos and antineutrinos. In contrast, in all the previous 
analyses the indication on the CP-phase 6 was indirect, since it arose from the combination of long-baseline (T2K) 
data and reactor data (Daya-Bay, RENO and Double-Chooz). Our analysis shows that the role of reactor data, albeit 
still significant, is now less pronounced in determining the estimate of the CP-phase 6 . This trend is likely to be 
reinforced in the future since the current LBL appearance data represent only a small fraction (~ 10%) of their total 
planned exposure. 

We recall that the LBL experiments are sensitive to CPV because they can probe the interference of two distinct 
oscillating frequencies. In the 3-flavor case, the two frequencies in question are related to the atmospheric squared- 
mass splitting and to the solar one. In the 4-flavor case, the interference occurs between the fast (averaged) oscillations 
induced by the new (9(eV^) splitting and the atmospheric one. As first shown in [SI], the new interference term has an 


I 


amplitude comparable to that of the standard 3-flavor one. Therefore, it is possible to extract from LBL experimenti 
some information on the enlarged CPV sector involved in the 3 - 1-1 scheme^ In particular, it was shown that the data 

preferred the common value dn ^ 614 - 7 r /2 (where di 3 s 6 ). Such a behavior was essentially driven by the slight 

mismatch existing between the two determinations of On achieved respectively with reactor and T2K data. 


We also recall that the LBL experiments are sensitive to the NMH because of the MSW effect 125 


2611 . which 


induce (tiny) modifications of the mass-mixing parameters in matter that depend on the sign of the atmospheric split¬ 
ting. In a sense, also this phenomenon may be regarded as a kind of interference, occurring between the MSW 
potential (whose sign is known) and the atmospheric squared-mass splitting (whose sign is unknown and to be de¬ 
termined). The same phenomenon is at the basis of the planned measurements of the NMH with atmospheric neutri- 


^For other recent studies addressing the issue of sterile neutrinos at LBL experiments see (HIIlIIl. Previous related works can be found 
in 

^In principle, also the atmospheric neutrinos may provide some information on the enlarged CPV sector, since they traverse very long distances 
in their path through the Earth. 
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3011 . where there is the advantage that the CPV interference effects are less important and there is no 


degeneracy among the determination of CPV and NMH, which instead notoriously plagues the LBL measurements. It 
is worth mentioning that the interference of the MSW matter potential with the atmospheric frequency is also present 
in the propagation of solar neutrinos (see for example j^), but unfortunately its amplitude is too small to extract any 
useful information from it (at least with the present data). 

Our 3-flavor fit evidences a consolidation of the previous hints, namely a slight preference for sin d < 0 and a mild 
preference (at more than 68 % C.L.) for the normal hierarchy. In a 3 - 1-1 framework, the data are able to constrain two 
(of the three) CP-phases (di 3 s d and 614 ), which exhibit a slight preference for the common value di 3 ^ 614 ^ -njl. 
Interestingly, in the enlarged four neutrino scheme the preference for the normal hierarchy found within the 3-flavor 
framework completely disappears. This occurs because the mismatch between the two estimates of 9\^ provided by 
reactor and LBL experiments (more pronounced in the IH case) can be “cured” with an appropriated choice of the new 
CP-phase ^ 14 . This circumstance indicates that light sterile neutrinos may constitute a potential source of fragility in 
the capability of the two LBL experiments of discriminating the neutrino mass hierarchy. 

The rest of the paper is organized as follows. In Sec. II we introduce the theoretical framework needed to dis¬ 
cuss the analytical behavior of the LBL and transition probabilities in the 3-flavor and 4-flavor 

frameworks. In Sec. Ill we present the results of the numerical analysis. Finally, in Sec. IV we draw our conclusions. 


2. Theoretical framework 

In the presence of a fourth sterile neutrino Vj, the flavor and the mass eigenstates are connected through a 4 x 4 
mixing matrix. A convenient parameterization of the mixing matrix is 


U - ^ 347 ^ 24 ^ 147 ^ 23 ^ 13^12 , 


( 1 ) 


where Rjj (Rij) represents a real (complex) 4x4 rotation in the (/, j) plane containing the 2 x 2 submatrix 
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in the (/, j) sub-block, with 


Cij — COS Gij 


Sij = sin Oij Sij = Sije 


( 2 ) 


( 3 ) 
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The parameterization in Eq. ([TJ enjoys the following properties; i) For vanishing mixing with the fourth state (614 = 
024 = 034 = 0) it returns the 3-flavor matrix in its usual parameterization, ii) The leftmost positioning of the matrix R ^4 


makes the vacuum ^ Ve conversion probability independent of 034 (and the related CP-phase ^ 34 ). iii) For small 
values of 0 i 3 and of the mixing angles involving the fourth mass eigenstate, one has |t/e 3 p - ij 3 , |t/e 4 p = ■S 14 (exact), 
\Un 4 \^ - *24 \ Ur 4 \^ ^ S 34 , with a clear physical interpretation of the new mixing angles. 

Fet us now come to the transition probability relevant for T2K and NOvA. From the discussion made in ISj], it 
emerges that the conversion probability can be approximated as the sum of three terms 


p4v pATM , pINT , pINT 


(4) 


The first term represents the positive definite atmospheric transition probability, the second term is related to the 
standard solar-atmospheric interference, while the third term is driven by the atmospheric-sterile interference. The 
probability depends on the three mixing angles ii 4 , S 24 - 0.15, which can be all assumed to be of the same order 
e and on the ratio of the solar and atmospheric squared-mass splittings a = ^ +0.03, which is of order e^. 

Keeping terms up to the third order, in vacuum, one finds 


)ATM 

“ 4^23^13 sin^ A, 

(5) 

)INT 

I 

^ 8ii3ii2Ci2i23C23(aA) sin Acos(A + ^ 13 ) , 

( 6 ) 

)INT 

II 

- 4 S 14 S 24 S 13 S 23 sinAsin(A + (5i3 - d^), 

(7) 


where A s Am^^LjAE is the atmospheric oscillating frequency, which depends on the baseline L and the neutrino 
energy E. 

The matter effects slightly modify the transition probability, introducing a dependency on the ratio 


_ y _ 2VE 

k AjMjj 

where 

V = ^/iGpNe 


(8) 


(9) 


is the (constant) matter potential along the neutrino trajectory. Both in T2K and in NOvA the value of v is relatively 
small, being v ~ 0.05 in T2K, and v ~ 0.17 in NOvA, where we have taken as a benchmark the peak energy of the two 
neutrino beams {E = 0.6 GeV in T2K, E - 2 GeV in NOvA). Therefore, it is natural to treat v as a small parameter 
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of order e. The probability in matter can be simply obtained (see the appendix in ySt] and the works 1132 . 
making in the vacuum expression the following substitutions 


33 


M)hy 


A^A™, (10) 

where and A'” represent the mixing angle S 13 and the atmospheric oscillation frequency in matter, which for small 
values of v, take the form 


^13 - (1+V)^13, 

A'" (l-v)A. 


( 11 ) 

( 12 ) 


It can be shown that at the peak energy, where |A| ^ njl, the substitutions in Eqs. (fTOl l - (fT^ produce third order 
corrections only in the atmospheric term, which reads 




lATM 


(13) 


while the same substitutions in the two interference terms produce only corrections of the fourth order. 
The swap of the NMH is parametrized by the replacements 


A —> -A, a -a, v -v. 


(14) 


One can observe that, while in general, the two interference terms are not invariant under a swap of the NMH, around 
the oscillation maximum |A| = 7 r/ 2 , such an invariance is approximately valid. Therefore, in practice the effect of 
the NMH swap is captured by the modification of the leading term in Eq. (flTt . Due to the change of sign of v, the 
transition probability (which acquires the dominant contribution from the atmospheric term) increases (decreases) 
with respect to the vacuum case in the NH (IH) case. NOvA is expected to be more sensitive than T2K to the NMH 
because of the larger value of the ratio v. 

Finally, we recall that the transition probability for antineutrinos is obtained from that of neutrinos by changing 
the sign of the MSW potential V and of all the CP-phases. This, for a given choice of the NMH, corresponds to the 
substitutions 


5i3 ^ -5i3, <5i4 ^ - 614 , V -V. 


(15) 
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In the NH case v > 0 for neutrinos and v < 0 for antineutrinos. According to Eq. (fOl) . in the NH case the leading 
contribution to the transition probability will increase (decrease) for neutrinos (antineutrinos). In the IH case the 
opposite conclusion holds. When passing from neutrinos to antineutrinos, the first interference term (around the 
oscillation maximum) changes sign, since it is essentially proportional to - sin ^ 13 . In contrast, the second interference 
term remains approximately invariant, since it depends on cos( 5 i 3 “ ^w)- 

Remarkably, for typical values of the mixing angles preferred by the current global 3+1 fits, the amplitude of 
the (atmospheric-sterile) interference term is almost identical to that of the standard (solar-atmospheric) interference 
term. As a consequence, a substantial impact on the regions reconstructed by the experiments T2K and NOvA in the 
pane of the two parameters [ 613 , di 3 ] is expected. In addition, a similar sensitivity to the two CP-phases di 3 and dn is 
expected in the combination of the two LBL experiments with the reactor data. 


3. Numerical Analysis 


In our numerical analysis we include the reactor experiments Daya-Bay and RENO and the two LBL experiments 
T2K and NOvA. The analysis of LBL data is slightly different for the two cases of three and four flavors, since in this 
last case there are appreciable oscillation effects not only at the far detector but also at the near detector, which have 
been taken into account as described in detail in Ist]. 

The analysis of the reactor experiments is performed using the total rate information and following the approach 


described in ll35n . Lor both experiments we have used the latest data yllj] based, respectively, on 621 live days 


(Daya Bay) and 800 live days (RENO). Since the electron neutrino survival probability probed by these experiments 
is independent of the CP-violating phases (standard and non-standard), their estimate of 0 i 3 is independent of them. 
We recall that such an estimate is extracted using the ratio of the event rates measured at the far and at the near sites. 
Since the fast oscillations induced by Aotj^ are averaged out at both detector sites, Daya Bay and RENO are not 
sensitive to 4-flavor effects. As a result, their estimate of ^13 is independent of the mixing angle 0 i 4 as long as it is 
allowed to vary in the range we are exploring. 

Concerning the LBL experiments, we use the T2K results of the ^ Vg appearance searches llj], which reported 
28 events with an estimated background of 4.92 events and the preliminary results of the ^ Vg searches presented 


at the EPS HEP conference 


ii3 


1 ih . which report 3 events with an estimated background of 1.68 events. Lor NOvA 


we use the preliminary results presented at Lermilab ||12|] . It must be noticed that the NOvA collaboration, prior 


to unbinding the data, decided to show the results obtained with two different event selection procedures (LID and 
LEM). The first is considered the primary method by the collaboration so we conservatively use the results obtained 
with such a method, which selects 6 events over an estimated background of 0.94 events. Lor completeness, we will 
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Figure 1: Left panels: regions allowed by the LBL experiments T2K and NOvA and by the 613 -sensitive reactor experiments for normal hierarchy 
(upper panel) and inverted hierarchy (lower panel). Right panels: regions allowed by their combination. The mixing angle 623 is marginalized 
away. The confidence levels refer to 1 d,o,f. (A^^ = 1.0,2.71). 

comment on the results that are obtained by using the secondary LEM selector, which identifies a considerably larger 
number (11) of events (the background being very similar to the LID case) in comparison to the primary LID method. 

In order to calculate the theoretical expectation for the total number of T2K Ve events and their binned spectrum 
in the reconstructed neutrino energy, we convolve the product of the v,j flux Q with horns operating in the -1-250 kA 
mode Oables provided on the T2K home page 1^ 1. the charged current quasi elastic (CCQE) cross-section (estimated 
from yj), and the —> Vg transition probability, with an appropriate energy resolution function. A similar procedure 
is followed for the antineutrino channel, where we have used the fluxes with horns operating in the -250 kA 


mode plotted in 11 111 . Concerning NOvA we used directly the product of the neutrino fluxes with the cross section 


(extrapolated at the far detector) provided in ||12|] . We have checked that our prediction for the binned spectra of events 


are in good agreement with those shown by the collaborations. We have performed the 3-flavor analysis both using 
the total rate information and the full spectrum, observing very small differences between the two methods. This is 
due to the limited statistics currently available, and to the effect of the smearing induced by the energy resolution. As 
explained in [St], in the 4-flavor case, from outside the collaborations, it is possible to consistently perform only a total 
rate analysis. Therefore, for homogeneity, also in the 3-flavor case, we report the results obtained with the total rate 
information. 


3.1. Results of the 3-flavor analysis 

In the 3-flavor analysis, the two mixing angles (0,3, 623 ) and the CP-phase (5,3 are treated as free parameters. 


taking into account the external prior sin^ 623 =0.51+ 0.05 provided by the disappearance measurement 1138 1 
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performed by T2K. For the atmospheric mass splitting we use the best fit value lAmjjl = 2.4 x 10“^ eV^ obtained in 
the same analysis. The solar mass-mixing parameters are fixed at the best fit value obtained in the global analysis |[^. 

Figure[T]shows the results of the analysis for the two cases of NH (upper panels) and IH (lower panels) in the plane 
spanned by the two variables [sin^ 20 i 3 , ^ 13 ], the atmospheric mixing angle 023 having been marginalized away. The 
left panels report the regions allowed by the combination of T2K (ve and Vg) and NOvA (vg) for the confidence levels 
68 % and 90% (1 d.o.f.). When possible [T2K (vg) and NOvA (Ve)], we have verified that the results of our analysis 
obtained for the single experiments in each channel return basically the contour plots presented by the collaborations. 
The narrow vertical band displayed in both panels represents the range allowed at 68 % C.L. for 0 i 3 by the reactor 
experiments. It is interesting to note that at the 68 % C.L., the contours determined by the combination of T2K and 
NOvA are closed regions around the best fit value bn ^ —7r/2. This means, that thanks to the information coming 
from the two channels (v^ ^ Vg and v^ —> Vg) the LBL experiments are starting to probe the CP symmetry in a direct 
way. This situation is qualitatively different from the preexisting one, where the extraction of the information on bn 
was indirect, i.e. not based upon a manifest observation of CPV. In fact, the previous hints in favor of bn - -7 t/2 
derived from the combination of T2K (vg) and reactor data. 

Fig. 1 shows appreciable differences between the two cases of NH and IH, which can be traced to the presence of 
the matter effects. As discussed in Sec. II, these tend to increase (decrease) the theoretically expected Vg rate in the case 
of NH (IH). The opposite is true for v^’s but their weight in the analysis is small, so the neutrino datasets dominate. 
In addition, as discussed in Sec. II, the NOvA Vg data are more sensitive than the T2K Vg data to the matter effects. 
Therefore, the neutrino datasets are far more relevant for what concern the sensitivity to NMH. More specifically, the 
following differences among the two hierarchies emerge, which will persist also in the 4-flavor analysis. The regions 
obtained for the case of IH: i) are shifted towards larger values of ^ 13 . ii) are slightly wider in the variable 0 i 3 with 
respect to those obtained in the NH case, and iii) tend to prefer (reject) values of sinbn < 0 (sinbn > 0) in a more 
pronounced way. 

The two right panels of Fig. 1 show that the combination of the reactor experiments with LBLs tends to further 
reinforce the preference for values of b ~ - 7 r/ 2 , disfavoring the case of no CPV (bi 3 = 0, tt) at roughly the 90% 
C.L. In addition, we note that the weak preference for the case of normal hierarchy tends to consolidate, being 
Tnh “Tih “ -1-3, to be compared with the result ~Xiu - “0-8 obtained in our previous work [SI]. However, the 
statistical significance of the indication is still small and below the 90% C.L. 

A final remark is in order concerning the two events selection methods used by NOvA. As already stressed above, 
we have conservatively adopted the results obtained with the primary selector, based on a likelihood identification 
(LID) method. By adopting the results obtained with the alternative method, based on a library event matching (LEM) 
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selector, which identifies a relatively larger number of Vg events (11 vs 6 ), we find that both the indication on the 
CP-phase ^13 and that on the NMH are slightly enhanced. More precisely, the rejection of the CP-conservation cases 
(di 3 = Q, 7 r) and the rejection of the inverted hierarchy both rise, respectively, at about the 2cr level and 1.5cr level. 
A combination of the results obtained with the two events selection methods is not trivial because they are strongly 
correlated. However, one can guess that such a combined analysis would provide results which are intermediate 
between those obtained using the two methods. 

3.2. Results of the 4-flavor analysis 

Figure |2] displays the results of the 4-flavor analysis for the case of NH. The four panels represent the regions 
allowed by T2K + NOvA in the usual plane [sin^ ^ 13 , ^ 13 ] for four different choices of the new CP-phase 614 . We have 
fixed the four-flavor parameters at the following values: s^^ = 5^4 - 0.025, sj^ - 0, ^34 = 0 and Ain^^ - 1 eV^. As 
a benchmark we also report the range allowed for On by reactors, which is identical to the standard case. A quick 
comparison of the four panels of Fig. |2] with the 3-flavor case (left upper panel of Fig.[TJ shows the noticeable impact 
of the 4-flavor effects on the allowed regions. The behavior of the curves can be qualitatively understood taking into 
account that the dominant contribution to the total rate comes from a region of the energies close to the first oscillation 
maximum, where A ~ njl, and that the neutrino datasets dominate over the antineutrino one. Inspection of Eq. dUl 
shows that the standard interference term is proportional to - sindi 3 . From Eq. O we see that for 614 = njl, the new 
interference term is proportional to sin ( 5 i 3 . Therefore, in this case the two terms are in opposition of phase and having 
similar amplitude tend to cancel out, making the wiggles structure almost to disappear (right upper panel). Vice 
versa, for 6\4 - -njl (right lower panel) the two inference terms have the same phase and the horizontal excursion 
of the “wiggles” is increased. In the two cases du - 0,n (left panels) the new interference term is proportional to 
+ cos ( 5 i 3 = + sin( 7 r /2 - fin) and thus it has a + 7 r /2 difference of phase with respect to the standard one. As a result, in 
such two cases, the behavior of the allowed regions is intermediate between the two cases 6\4 - (- 7 r/ 2 , 7 r/ 2 ). 

Eigure[3]shows the allowed regions for the case of IH. Eor each choice of 614 , when going from NH to IH one finds 
similar changes to those already discussed for the 3-flavor case (shift of On towards larger values, a wider allowed 
range for On and a more pronounced preference for sindi 3 < 0). The changes induced by the hierarchy swap are 
similar to the 3-flavor case because the matter effects enter in the same way in the 3-flavor case and in the 4-flavor 
one, i.e. modifying the dominant atmospheric term in the transition probability [see Eq. mi 

It is interesting to note how, in the presence of 4v effects, a better agreement among the two estimates of On 
derived from reactors and LBLs can be obtained. In particular, this occurs for du ^ -njl. As we have discussed for 
the 3-flavor case, the mismatch of the 613 estimates from LBLs and reactors tends to disfavor the inverted hierarchy. 
The same is no more true in the 4-flavor scheme, since the two estimates can be brought in agreement for an appro- 
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Figure 2: Regions allowed by the combination of T2K and NOvA for four representative values of the CP-phase 614 . Normal hierarchy is assumed. 
The mixing angle 623 is marginalized away. The vertical band represents the region allowed by reactor experiments. The confidence levels are as 
in Fig,[T] 



Figure 3: Regions allowed by the combination of T2K and NOvA for four representative values of the CP-phase ^ 14 . Inverted hierarchy is assumed. 
The mixing angle 623, is marginalized away. The vertical band represents the region allowed by reactor experiments. The confidence levels are as 
in Fig.^ 


priate choice of the new CP-phase (^14 ^ -7r/2) (see the right bottom panel ofO. This circumstance indicates that 
light sterile neutrinos may constitute a potential source of fragility in the capability of the two LBL experiments of 
discriminating the neutrino mass hierarchy. 

As a last step in our 4-flavor analysis, we perform the combination of LBLs with reactors. In this more general 
analysis, we treat the two mixing angles ( 0 i 3 , 623 ) and the two CP-phases (^ 13 ,^ 14 ) as free parameters, while fixing the 
remaining 4-flavor parameters at the same values used before; s\^ — s\^ - 0.025, S 34 = 0, ^34 = 0 and AiMj^ = 1 eV^. 
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Figure 4: Regions allowed by the combination of the LBL experiments and the reactor experiments for the case of normal hierarchy. The mixing 
angle 023 is marginalized away. 


We have checked that the impact of non-zero 634 (and consequently of the associated CP-phase ^ 34 ) is negligible when 
considering values of 634 below the curi'ent upper bounds. 

Similarly to the 3-flavor case, in the LBL - 1 - Reactor combination, the (CP-phases independent) estimate of 6*13 
provided by the reactor experiments selects those subregions of the LBL bands which have a statistically significant 
overlap with such an estimate. These, in turn, correspond to allowed regions in the plane [( 5 i 3 , 5 i 4 ] spanned by the 
two CP-phases. Figures|4]and|5]display such regions for the two cases of NH and IH, together with the 2-dimensional 
projections having as one of the two variables the mixing angle 613 . As expected, there is a similar sensitivity to 

both CP-phases with a preference for values of di 3 ~ di 4 - njl. The case dn = 0 is disfavored at a slightly lower 

confidence level in comparison with the 3-flavor case. This is imputable to the degeneracy among the two CP-phases. 


4. Conclusions and Outlook 


We have investigated the impact of the latest data released by the two long-baseline experiments T2K and NOvA. 
Specifically, we have included in our analysis the first Vg appearance data from T2K and the first Ve appearance results 
from NOvA. We have considered the implications of such new results both for the standard 3-flavor framework and 
for the non-standard 3 h- 1 scheme involving one sterile neutrino species, which is the most popular scenario invoked 


to explain the short-baseline anomalies (see 


mmy 


The 3-flavor analysis shows a consolidation of the preexisting trends, namely a slight preference for sin 5 i 3 < 0 
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Figure 5: Regions allowed by the combination of the LBL experiments and the reactor experiments for the case of inverted hierarchy. The mixing 
angle ^23 is marginalized away. 


disfavoring the cases of CP conservation (5 = O.tt) at a confidence level close to 90% C.L., and a mild preference 
(at more than 68 % C.L.) for the normal hierarchy. These hints are slightly reinforced if one adopts the NOvA results 
obtained with the secondary event selection method (LEM). In a 3+1 framework, two CP-phases (( 5 i 3 s 6 and ( 5 i 4 ) 
are constrained, which exhibit a preference for the common value ^13 ^ 614 ^ - 7 r/ 2 . Interestingly, in the enlarged 
four neutrino scheme the preference for the NH found within the 3-flavor framework completely disappears because 
the mismatch between the two estimates of On from LBL’s and reactors can be lifted (an intriguing circumstance) if 
the new CP-phase du lies around its best value -n/l. This fact indicates that light sterile neutrinos may constitute a 
potential source of fragility in the capability of the two LBL experiments of discriminating between the two neutrino 
mass hierarchies. 

Our results show once again that LBL experiments can give an important contribution in the context of sterile 
neutrino searches, which should not considered only a short-baseline “affair”. In particular, no information on the new 
CP-phase ^14 can be extracted from SBL setups, where the interference between the new splitting and the atmospheric 
one cannot develop (or is completely negligible). Therefore, the two classes of experiments (SBL and LBL) are 
complementary and synergic in constraining the 4-flavor parameter space. While a discovery of a sterile neutrino can 
come only from the observation of the oscillating pattern (in energy or/and space) at the new SBL experiments, the 
full exploration of the 4v model will become possible only with the contribution of LBL setups, which are sensible 
interferometers able to probe the new enlarged CPV sector. Also the (present and future) atmospheric neutrino data 
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may help to shed light on this issue, since these neutrinos traverse very long baselines and may present some sensitivity 
to the new CPV phenomena. 

Finally, we would like to underline the other side of the coin, namely the fact that the presence of sterile neutrinos 
tend to decrease the robustness of the conclusions reached in the simple 3-flavor framework. The present data already 
show that this is the case for the current hint in favor of NH, which disappears in the 3-1-1 scheme. In this respect, 
it would be very interesting to study how the NMH discrimination capabilities of future data expected to come from 
T2K and NOvA, and from the planned LBL experiments (DUNE, LBNO and T2HK) will be affected by the inclusion 
of the sterile neutrino effects. 
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